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TRANSIENT  RESPONSE  OF  LONG  STRUCTURAL  MEMBERS 


By  David  Hsu^  and  Walter  Pilkey^,  M.  ASCE 


INTRODUCTION 

The  dynamic  response  of  long  structural  members  has  usually  been 
approached  using  wave  propagation  methodology.  JJith  such  methods  an  arbitrary 

structural  member,  with  variations  in  geometry  and  material  and  complicated 

> 

loading,  is  a  challenge  to  analyze  for  transient  response.  A  simple  pro¬ 
cedure  is  proposed  here  which  reduces  the  long  member  analysis  to  that  of 
a  finite  length  member  analysis.  This  permits  a  long  member  with  arbitrary 
loading,  material,  and  geometry  to  be  studied  using  the  integration  techniques 
commonly  employed  for  finite  members. 

The  proposed  procedure  first  requires  the  reduction  of  the  governing 
partial  differential  equation  to  an  ordinary  differential  equation  in  the 
spatial  variable.  The  time  variation  can  be  removed  using  a  finite  difference 
discretization  or  with  a  transform.  Here,  we  employ  a  Laplace  transform.  The 
resulting  ordinary  differential  equation  in  the  spatial  variable  can  be  solved 
for  a  long  member  using  the  transfer  matrix  procedure  for  infinite  beams  of 

Ref.  (1).  This  method  will  accept  members  with  arbitrary  loading 
and  changes  in  cross-sectional  material  and  geometric  properties.  The 

final  solution  in  the  time  domain  is  then  obtained  by  inverting  the 
solution  in  the  transform  domain.  Usually  the  inverse  transform  must  be 
computed  numerically. 
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BASIC  THEORY 


The  procedure  will  first  be  illustrated  using  the  longitudinal  motion 
of  a  bar.  The  governing  equations  of  motion  for  the  axial  extension  of  a 
bar  are  (2) 


3u 

3x 
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+  k  u 
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(1) 


in  which  u  ■  Axial  displacement,  P  ■  axial  force,  E  -  modulus  of  elasticity, 

A  *  cross-sectional  area,  p  *  mass/length,  px  ■  distributed  axial  force, 
o  *  coefficient  of  thermal  expansion,  AT  *  change  of  temperature,  kx  «  elastic 
foundation  modulus.  To  simplify  the  equations,  ignore  the  aAT,  and  kxu  terms. 
Then  Eq3.  (1)  can  be  combined  to  give  the  second  order  equation 

3^u  1  3^u  Px(x,C) 

F”c2F*  EA -  (2) 


where  c  ■  yEA/ p  , 

Consider  a  very  long  bar,  such  as  a  semi- inf inite  bar  (x>0) .  We  shall  assume  that 
Px(x,t)  •  0  for  x  >  a  and  that  all  loadings,  changes  of  materials,  variations 
of  cross-sectional  area,  etc.  occur  only  to  the  left  of  x  ■  a.  In  other  words, 
we  are  treating  Eq.  (2)  in  two  parts: 
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where  0  ■  a  <  a,  <  a0  <...<  a  "a 
oiz  n 

and 
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To  solve  Eqs.  (3)  and  (4)  with  initial  conditions  and  boundary  conditions. 


A 


use  the  Laplace  transformation  on  the  independent  variable  t,  so  that  the 
partial  differential  equations  are  reduced  to  the  ordinary  differential 
equations 
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-  -iy  s2u  -  su(x,  0)  -  u(x,  0)  - 
dxz 


rix 

rr~  a<  1  x  i  aa  7  i  -  i,. . .  ,n  (5) 


— 2  s2u  -  su(x,0)  -  u(x,0) 


x  >  a 


where  u  »  Laplace  transform  of  u,  P  «  Laplace  transform  of  P.  The  left-end 
boundary  condition  at  x  *  0  is  given,  while  the  boundary  condition  at  the 
infinite  end  requires  a  finite  value  of  u  for  x  *  =». 

For  u(x,0)  *  u(x,0)  *  0,  Eq.  (6)  has  the  solution 
-1  (x-a)  1  (x-a) 

u(x,s)  »  ae  +  8e  (7) 


and  the  bounded  response  requirement  at  infinity  demands  that  3  ■  0.  Define 

—  —  T  T 

the  state  vector  z(x)  *  [u  P  1]  .  Also,  define  C  *  [a  0  1]  .  For  x  >  a,  note 

that 


To  compute  C  in  terms  of  the  state  vector  at  a,  substitute  x  ■  a  in  Eq.  (8) 
and  form 
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For  x  ^  a,  use  the  usual  transfer  matrix  method 

z(x)  -  Ui(x,ai_1>  •••  2(°)  ai  >  x  >  ai_1  (10) 

with  U.  being  the  transfer  matrices  in  [a^^.a^]  obtained  from  Eqs.  (5).  Use 
of  the  transfer  matrix  method  permits  the  inclusion  of  arbitrary  loadings  and 
variations  in  material  and  geometry. 

One  of  the  state  variables  at  x»0  (u(0,s)  or  P(0,s))  is  known  from  the 
boundary  condition  at  x*0.  The  other  state  variable  at  x*Q  can  be  determined 
by  equating  z(a)  ■  uQ(a»an_i^  Un-l^an-l,an-2^  *"  (0)  -  U(a,0)z(0) 

and  z(a)  ■  Q(a,a)C,  i.e, 

U(a,0)z(0)  *  Q(a,a)C  (11) 

This  expression  fixes  a  and  then  is  used  to  evaluate  the  unknown  state 
variable  at  x  »  0.  Now  that  both  u(0,s)  and  P(0,s)  are  known,  i.e.,  z(0) 
has  been  found,  the  state  vector  z(x)  can  be  calculated  using 

for  ajL  _>  x  _>  a^  i  -  1, 2, . .  .  ,n 

z(x)  -  U1(x,ai-1)Ui^1(a1_1,ai_2)  ...  U1(a1,0)z(0)  (12) 

for  x  ^  a 

z ( x)  «  Q(x,a)C  -  Q(x,a)Q’1(a,a)z(a)  -  Q(x,a)Q  1( a,a)U(a, 0) z(0)  (13) 

The  solution  u(x, s)  is  then  transformed  back  to  the  time  domain  through  the 
formula 

u(x,t)  -  /*u(x,s)estds 
*T 

If  the  integral  of  Eq.  (14)  is  too  difficult  to  evaluate  analytically, 
the  discrete  fast  Fourier  transform  routine  (3)  can  be  used  for  a  numerical. 


approximate  evaluation. 
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If  the  bar  is  semi- inf  inite  to  the  left  (x  0) ,  o  rather  than  8  should 
be  set  to  0,  and  Eq.(ll)  would  be  replaced  by 

U(b,0)z(0)  -  Q(b,b)D  with  r  *  '0  S  1]T  (15) 

Here  b  is  the  point  to  the  left  of  which  no  variation  in  material,  geometry, 
or  loading  can  occur.  For  bars  which  are  infinite  at  both  ends,  Eq.  (15) 
is  changed  to 

Q(a,a)C  -  U(a,b)Q(b,b)D  (16) 

from  which  we  determine  a  in  C  and  0  in  D,  or,  equivalently,  the  state  vector 
at  x*0. 

BEAMS 

For  the  dynamic  response  of  beams,  the  4th  order  governing  homogeneous 
equation  after  the  Laplace  transformation'  is 

fyr  +  U  -  n)  0  +  (  -  Cn)  S  -  0  (17) 

where  5  »  (P  -  pr2s2  -  k*)/El,  n  »  (k  +  ps2)/GAs,  X  *  (k  +  ps2)/EI 

The  initial  conditions  are  taken  to  be  zero  in  Eq.  (17).  The  nomenclature  of 
Ref.  (2)  has  been  used.  The  quantities  Z,  n,  and  X  are  complex  because  of  the 
presence  of  s.  Let  u(x,s)  ■  e  rx  and  substitute  it  into  Eq.  (17). 

r4  +  (z  -  r^r2  +  (X  -  ?n)  -  0 
Since  the  coefficients  are  complex,  denote 

p  +  qi  -  U  -  n)2  -  4<x  -  cn),  a  +  bi  -  (5  -  n) 


(18) 


(19) 


Hence  the  general  solution  for  Eq.  (17)  is 


u(x,s) 
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Since  h^  >0,  h2  >  0  for  any  value  of  s,  and  C2  must  be  set  to 
zero  in  order  for  u  to  be  bounded  at  x  «  •.  It  is  readily  shown  that  for 
the  state  vector  z(x)  *  [u  9  M  V]  and  C  ■  [C^  C2  C_  C^] 


The  two  unknown  state  variables  at  x*0  can  be  determined  from  U(a,0)z(0)  ■  Q(a,a)C. 

The  rest  of  the  procedure  remains  the  same  as  that  for  the  axial  motion  of 

bars. 


EXAMPLE 

As  an  example  of  this  technique,  the  stepped  semi-infinite  bar  of  Fig. 

1  was  analyzed.  The  applied  loading  is  P(t)  -  *>0te"’t-  7116  initial  conditions 

are  u(x,0)  «  0}  u(x,0)  »  0,  the  boundary  condition  at  the  left  end  is 
P(0,t)  *  0,  and  u(a>,t)  is  bounded. 
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The  transfer  matrix  U^x.a^,  for  a2>.  xj^,  is  obtained  from  Eq.  (22)  by  re¬ 
placing  x,  Ex,  ^  by  x-a1,  c 2,  E 2,  A^  respectively.  The  quantity  Q(x,a2), 
xZ.a2»  ta^en  from  Q(x,a)  by  replacing  a,  c,  E,  A  by  a2,  Cy  Ey  A^,  respectively. 
From  U(a2,0)z(0)  -  Q(a2,a2)C,  the  initial  parameter  u(0,s)  is  found  to  be 


u(o,s)  - j- 


E_A_c- 

P_/a  —  ~  sinh|^(a2  -  a^  +  cosh  ~  al) 


V1 


E2A2C3 


H 


(24) 


with 


a  sE*  A-  s 

h  H  -  — - sinh^  (a2  -  a^  +  — - —  cosh  -^(a2  “  al)sinh  ^ai 


'2 
‘VS 


sc2E.A2E-A- 

+  — :  '  cosh  f-  (a,  *  a.  )  .  cosh  7  a..  +  •  ■  v  t -  sinh-  (a2  -  a^ 

c2  2  1  C1  1  C1C3E2A2  2 


8 


For  a  more  conplicated  problem  the  evaluation  of  u(0,s)  and  u(x,s)  would  be 
handled  numerically.  The  transfer  matrix  method  can  now  be  used  to  provide 
u(x,s)  for  any  x. 

Figure  1  shows  the  displacement  at  x  *  15  for  the  bar  with  A^  *  1,  A2  •  2, 

A3  *  3»  *  0*002,  p2  -  0.004,  p3  -  0.006,  Ej_  -  E2  =*  Ej  -  107,  Pq  -  20,  ^  -  10, 

a2  »  20. 

SUMMARY  AND  CONCLUSIONS 

A  technique  is  sec  forth  by  which  the  dynamic  response  can  be  found 
for  long  structural  members  with  arbitrary  axial  variations  in  applied  loading, 
geometry,  and  material.  This  simple  procedure  reduces  the  long  member  analysis 
to  that  of  a  member  of  finite  length. 
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